typeset using JPSJ.sty <ver.l.Ob> 



Critical Phenomena of Ferromagnetic Transition in 
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Critical phenomena of ferromagnetic transition at finite temperatures are studied in double- 
exchange systems. In order to investigate strong interplay between charge and spin degrees of 
freedom, Monte Carlo technique is applied to include fluctuations in a controlled and unbiased 
manner. By using finite-size scaling analysis, critical exponents and transition temperature are 
estimated for a model with Ising spin symmetry in two dimensions. The obtained exponents are 
far distinct from the mean-field values, but consistent with those of spin models with short-range 
exchange interactions. The universality class of this transition belongs to that of short-range 
interaction with the same spin symmetry. We also discuss the case for three dimensions. The 
results are compared with experimental results in perovskite manganites which show colossal 
magnetoresistance. 
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Perovskite manganites which show the colossal mag- 
netoresistance have revived interests in the double- 
exchange (DE) mechanism.El^ Recent intensive experi- 
ments have improved the precision of data remarkably,© 
which now makes an issue of detailed comparisons be- 
tween experiment and theory in a quantitative way. Al- 
though theoretical survey on the DE systems has a long 
history, many aspects at finite, temperatures have not 
been fully understood thus farJ3' For instance, the criti- 
cal temperature of the metallic-ferromagnetic transition 
has been controversial until quite recently.o' 

What makes it difficult to describe the system is strong 
spin fluctuation. The system at issue is a strongly- 
correlated electron system since Hund's-rule coupling be- 
tween conduction electrons and localized spins is much 
larger than the kinetic bandwidth. Electrons move 
around and align localized spins in a parallel way to 
gain the kijaetic energy, which is an origin of the DE 
interaction.liP At finite temperatures, especially near the 
transition, spin fluctuations strongly modify the kinetics 
of electrons and vice versa through the strong interplay 
between charge and spin degrees of freedom. 

One of the most fundamental and challenging prob- 
lems is critical phenomenon of the ferromagnetic tran- 
sition in this system. The DE interaction has a distin- 
guishable property compared to the ordinary exchange 
interactions in localized-spin models such as the Heisen- 
berg model. Since the interaction is originated from the 
motion of conduction electrons, the effective range of 
the interaction may depend on kinetics of electrons and 
hence on temperatures crucially through the strong fluc- 
tuations. It is not clear how the interaction is renormal- 
ized near the critical point and what universality class 
governs the magnetic transition. 

To be more specific, we consider a situation where spin 



configuration forms ferromagnetic domains as illustrated 
in Fig. 1, as an example. In the limit of strong Hund's- 
rule coupling, conduction electrons are completely par- 
allel to localized spins and confined within the domains. 
In this system, the total energy is given by a sum of 
energy eigenvalues e^, i^tot = '^Su, where summation 
is restricted to occupied states. Energy eigenvalues e^, 
for electrons confined within domains strongly depend 
on the sizes as well as the shapes of the domains. In or- 
der to describe i?tot as a function of spin configuration, 
it is necessary to introduce long-range two-spin interac- 
tions as well as multiple-spin interactions which depend 
on sizes and shapes of ferromagnetic clusters. This is a 
marked contrast to the ordinary spin models with short- 
range exchange interactions, where the total energy de- 
pends only on surface volumes of domain boundaries. 
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Fig. 1. Schematic picture of the double-exchange model in the 
limit of strong Hund's-rule coupling. White and gray areas 
describe ferromagnetic domains where electrons are localized 
within. 

The size of domains grows up and the shape fluctuates 
strongly as the system approaches to the critical point. 
In the DE systems, it is difficult to describe the renor- 
malization of the effective magnetic interaction near the 
transition, whether it is renormalized to a short-range 
interaction, or the long-ranged and multiple-spin interac- 
tion become relevant to cause mean-field-like transition 
through suppression of fluctuations. Thus, it is highly 
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nontrivial and interesting to determine the universahty 
class of the ferromagnetic transition of the DE origin. 

Recently, we studied the critical phenomena in the 
three-dimensional (3D) DE modelJj-' The critical expo- 
nents are examined by using the finite-size scaling anal- 
ysis on the Monte Carlo (MC) data which is in principle 
an unbiased method. Our data scale more consistently 
to the short-range Heisenberg exponents than to those 
of mean-field results. However, the data are not pre- 
cise enough to determine the universality class in detail, 
mainly because the system size is still not large enough 
due to the limitation in computational time. 

In this letter, we examine critical phenomena of the 
metallic ferromagnetism in the DE systems, through a 
precise estimation for critical exponents. For this pur- 
pose, we study a simplified DE model which still cap- 
tures the essence of the DE mechanism. By analyzing 
the MC data using the finite-size scaling, we estimate 
the critical exponents as well as the critical tempera- 
ture. We compare the results with those for spin models 
with short-range type interactions as well as mean-field 
values. 

In order to perform calculations on a larger scale more 
precisely, we study minimum models which include the 
DE mechanism and show the ferromagnetic transition at 
finite temperatures. One of such models is the DE model 
where electrons are coupled to spins with Ising symmetry 
in two dimensions. The Hamiltonian is given by 



<JJ>,<T 



JhY^ctI-S,, (1) 



where Cia (c|^) annihilates (creates) a a-spin electron 
at site z, Jh is the Hund's-rule coupling, is the z- 
component of the spins of conduction electrons, and S 
describes the localized Ising spin which takes = ±1. 
The summation in the first term is taken for nearest- 
neighbor sites. For simplicity, we consider the limit of 
strong Hund's-rule coupling, Jh oo. Then the elec- 
tron hopping between sites i and j is allowed only when 
the localized spins are parallel, that is. Si — Sj. The 
model becomes equivalent to a spinless-fermion model 
with hopping integrals which depend on the spin config- 
urations: 



1 -|- SiSj) c\c- 



(2) 



where Ci and c| are the spinless-fermion operators. Note 
that model (||) is simplified compared to the original 
model. Nevertheless, it still includes the essence of the 
DE mechanism; the kinetics of electrons is strongly cor- 
related with spin degrees of freedom and energetically 
favors ferromagnetism. p. 

We study model (^ by an improved MC methodJj^ 
The most time-consuming part in the standard MC algo- 
rithm is the exact diagonalization of the fcrmion Hamil- 
tonian for each spin configuration.^' In the present al- 
gorithm, it is replaced by moment expansion of density 
of states. This reduces the cost of computational time 
considerably, while the truncation error falls off exponen- 
tially. The algorithm has another advantage to reduce 



computational time; the moment expansion is performed 
on parallel computers efficiently. 

In the following, the ferromagnetic transition in model 
(|^) is studied at the electron filling x = {J^i cjci) = 0.5. 
Here, the bracket denotes the thermal average for the 
grand canonical ensemble. A closed-shell condition is 
chosen to obtain rapid convergence of MG samplings and 
systematic extrapolation of system sizes.© We study the 
system with square lattice up to 20 x 20 under the pe- 
riodic boundary condition in one direction and the an- 
tiperiodic boundary condition in the other direction. We 
take the half-bandwidth W of noninteracting electrons as 
an energy unit; = 4t in two dimensions. The moment 
expansion of density of states is taken up to 20th order. 
We confirm that the truncation error is negligible small 
in the region T/W > 0.05 in model (||) in comparison 
with the results by exact diagonalization method. We 
have typically run 100,000 MC samplings for measure- 
ments after 10,000 MC steps for thermalization. Mea- 
surement is divided into five bins to estimate the statis- 
tical error by the variance among the bins. 

In order to determine the critical exponents and the 
critical temperature Tc, we use the finite-size scaling 
analysis on the ferromagnetic component of the spin 
structure factor. The spin structure factor is defined by 

Sik)^l^Y.(^,S,)e'''-^^', (3) 

where N is the total number of sites. If we assume the 
hyperscaling hypothesis, the ferromagnetic component 
S{k — 0) satisfies the scaling relation near the ferromag- 
netic transition as 



S'(0)L''-2 = fieL^/"), 



(4) 



where e = (T — Tc)/Tc, L is a linear dimension of the sys- 
tem and / is the scaling function. Here rj and v are the 
critical exponents, by which we can determine the uni- 
versality class under the hyperscaling hypothesis. For 
given exponents and Tc, we plot S{0)L'^~^ as a function 
of eL^/'^ and fit the data by a polynomial function /gt 
which approximately gives the scaling function /. The 
best estimate for exponents and is given by minimiz- 
ing the systematic deviations of the data from the scaling 
function. We define the deviation function by differences 
between the fit and the MC data as 



(5) 



where ?/„ = Sn{Q)L^'^ and Xn = ^uLI!'' . Here the sub- 
script n labels MC data for different parameters T and 
L. We apply the scaling of eq. (^) on MC data for the 
systems with L = 12, 14, 16, 18 and 20 for 0.05 < T/W < 
0.065. We use the fitting function /fit of the polynomials 
up to 12th order. For convenience of later discussions, 
hereafter we use the exponent (3 which is related with 77 
and vhy (3 — riv/2 in two dimensions. The best estimate 
which minimizes R is obtained at /3 = 0.09, v = 0.9 and 
Tc/W = 0.058. The scaling behavior at these values is 
shown in Fig. 2. All the MC data are scaled well to a 
single universal function within the errors. 
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Fig. 2. The best-fit result of the scaUng plot for the Monte Carlo 
data. The gray curve in the figure is the polynomial fit. 

Errors of the critical exponents are estimated in the 
following way. Since the MC data have statistical errors, 
the deviation function (|^) also has an error 6R. In the 
parameter space of iy,Tc), there exists a region where 
the values R are not distinct from its minimum value 
within the statistical errors, by which we define the re- 
gion for the estimate of {f3, v, T^). Namely, we determine 
the error region by the condition 



R-SR<R* + SR* 



(6) 



where R* and 6R* are the minimum value of R and its 
error at the best-fit point, respectively. The gray regions 
in Fig. 3, including hatched areas, show the results by 
making projections of the error region onto /3-i^ and P-Tc 
planes. 
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Fig. 3. Estimated values of (a) the critical exponents u and /3; 
and (b) the critical temperature Tc and /3. The circle, square 
and diamond in the figures denote the best-fit result, the Ising 
exponents and the mean-field exponents, respectively. The gray 
areas, both with and without hatches, exhibit the estimates by 
the scaling analysis of the spin structure factor. The hatched 
areas are the results by using the estimate for Tc by the Binder 
parameter in addition. See text for details. 

As an alternative method to estimate Tc indepen- 
dently, we perform the Binder parameter analysis.™ The 
Binder parameter is defined by g = 1 — (9^) /?>{q^)'^ where 
the order parameter q is taken as Si in this case. 
For T > Tc {T < Tc), the Binder parameter decreases 
(increases) as the system size N increases. Thus Tc is 



determined by a crossing point of g for various N. 

Figure 4 shows the temperature dependence of g for 
different system sizes. The inset displays the system-size 
dependence. The Binder parameter seems to cross at 
T/W ~ 0.058. The data at T/W = 0.058 are almost 
independent of N as shown in the inset. We estimate 
Tc/W = 0.058±0.001. This is consistent with the results 
in Fig. 3. Using this value of Tc as an additional criterion 
in the analysis of eq. we obtain the estimates for the 
critical exponents as shown by the hatched regions in 
Fig. 3. As a conclusion from combined analysis of finite- 
size scaling and Binder plots, we obtain f3 — 0.09 ±0.08, 
v = 0.9± 0.3 and Tc = 0.058 ± 0.001. In Table I, we 
show a summary of the exponents. 
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Fig. 4. Temperature dependence of the Binder parameter. The 
inset is system-size dependence; the data are for T/W = 
0.061,0.06,0.059,0.058,0.057,0.056,0.055 and 0.054 from bot- 
tom to top. The curves are guides to the eye. 



Table I. Comparison of the critical exponents. 

/3 



present model 
2D Ising model 
mean-field 



0.09 ± 0.08 
0.125 
0.5 



0.9 ±0.3 
1 

0.5 



The exponents and Tc in Fig. 3 are crosschecked by 
the scaling behavior of the magnetization. The magne- 
tization is calculated by m = [S{k = 0)/N]^^^. The 
temperature dependence of the magnetization is shown 
in Fig. 5. The values in the thermodynamic limit are 
obtained by system-size extrapolation of the spin struc- 
ture factor, which is shown in the inset. MC data are 
scaled well to which is expected in the ordered 

phase in the limit of large N. We fit the magnetiza- 
tion by assuming the scaling relation m oc (Tc — T)^ 
and obtain Tc/W = 0.0584 ± 0.0005 and the exponent 
(3 = 0.11 ± 0.01. If we assume that these data are within 
the critical region, these values of Tc and /3 are all con- 
sistent with the results in Fig. 3. Note that in this tran- 
sition it is difficult to determine Tc and /3 by the magne- 
tization data alone because of the small value of (3. 

Now, we discuss the results in Fig. 3. For compar- 
ison, we plot the critical exponents of the Ising model 
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with nearest-neighbor interaction in two dimensions; /3 = 
0.125 and v = lBa> These values are included in our es- 
timated regions and very close to the best-fit result. On 
the contrary, the critical exponents of the mean-field ap- 
proximation, f3 — — 1/2, are far outside of the area 
as shown in the figure. These results indicate that the 
universality class of the ferromagnetic transition in our 
simplified model (^) is the same as that in models which 
have short-range two-spin interactions but clearly dis- 
tinct from that of the mean-field approximation. 
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Fig. 5. Temperature dependence of the magnetization. The 
curve is the least-square fit to m oc {T — Tc)^ ■ The in- 
set shows system-size extrapolation of the spin-structure fac- 
tor. The lines are the least-square fits to 1/7V. The data are 
for T/W = 0.05, 0.051, 0.052, 0.053, 0.054, 0.055, 0.056 and 0.057 
from top to bottom. 

Let US now consider the 3D DE model. The nature of 
the magnetic interaction which is mediated by itinerant 
electrons are common. Estimates for the critical expo- 
nents in the 3D DE model are similar to the shoitt-range 
Heisenberg type rather than the mean-field onea^ Com- 
parison between the present result and those for 3D DE 
model strongly suggests the following: The long-range 
and multiple-spin parts of the magnetic interactions in 
the DE models seem to be renormalized to be irrelevant. 
The universality class of the ferromagnetic transition in 
the DE systems should belong to that of the short-range 
interaction with the same spin symmetry. 

Concerning experiments in manganites, there have 
been many studies to estimate the critical exponents 
of the ferromagnetism. Up to now, however, the ex- 
perimental estimates for exponents are still controver- 
sial, irply ding—I. hose for short-range Heiseiiberg-.inter- 
action,BBByM the mean-field values,eM'tZ|) and 
those which camiat be classified into any universality 
class ever known.tS' 

From our numerical results, it is suggested that the ex- 
ponents for the short-range Heisenberg universality class 
should be obtained when the transition in real materials 
is ascribed to the DE mechanism as a major origin. Fur- 
ther experimental studies are desired to clarify the criti- 
cal exponents of the real materials, and examine whether 
the transition follows the single-parameter scaling in the 
universality class of the short-range type. These include 



sample-quality refinements, especially in homogeneities 
of Tc, as well as appropriate treatments for the fitting 
within critical regions. The experimental results should 
cast a crucial test whether the DE mechanism alone can 
explain experimental results or not. 

In summary, we have studied critical phenomena of 
the ferromagnetic transition in the simplified double- 
exchange system numerically. The model which has Ising 
spin symmetry has been studied in the limit of strong 
coupling on square lattices. We have applied an im- 
proved Monte Carlo algorithm which reduces computa- 
tional time considerably. Critical exponents and critical 
temperature are estimated by the finite-size scaling anal- 
ysis on the Monte Carlo data. The critical exponents 
P = 0.09 ± 0.08 and v = 0.9 ± 0.3 are consistent with 
those of models with short-range interactions, but dis- 
tinct from those of the mean-field approximation. The 
universality class of this transition belongs to the short- 
range type with the same spin symmetry in the model. 

The authors thank H. Nakata for helpful support in de- 
veloping parallel-processing systems. The computations 
have been performed mainly using the facilities in the 



AOYAMA-I- project ( http:/ /www. phys.aoyama.ac.jp/ 



"aoyama-l-) and in the Supercomputer Center, Institute 
for Solid State Physics, University of Tokyo. This work 
is supported by "a Grant-in-Aid from the Ministry of 
Education, Culture, Sports, Science and Technology" . 



[1] C. Zener: Phys. Rev. 82 (1951) 403. 

[2] A. P. Ramirez: J. Phys.: Condons. Matter 9 (1997) 8171, and 
references therein. 

[3] N. Furukawa: in Physics of Manganites, eds. T. Kaplan and 
S. Mahanti (Plenum Publishing, New York, 1999), and refer- 
ences therein. 

[4] Y. Motome and N. Furukawa: J. Phys. Soc. Jpn. 69 (2000) 
3785. 

[5] Y. Motome and N. Furukawa: J. Phys. Soc. Jpn. 68 (1999) 
3853. 

[6] S. Yunoki, J. Hu, A. Malvezzi, A. Moreo, N. Furukawa and 

E. Dagotto: Phys. Rev. Lett. 80 (1998) 845. 
[7] K. Binder: Z. Phys. B 43 (1981) 119. 

[8] B. Kaufmann and L. Onsager: Phys. Rev. 76 (1944) 1244. 
[9] J. Stephenson: J. Math. Phys. 5 (1964) 1009. 
[10] R. H. Heffner, L. P. Le, M. F. Hundley J. J. Neumeier, G. 

M. Luke, K. Kojima, B. Naomichi, Y. J. Uemura, D. E. 

MacLaughlin and S-W. Cheong: Phys. Rev. Lett. 77 (1996) 

1869. 

[11] M. C. Martin, G. Shirane, Y. Endoh, K. Hirota, Y. Moritomo 

and Y. Tokura: Phys. Rev. B 53 (1996) 14285. 
[12] S. E. Lofland, S. M. Bhagat, K. Ghosh, R. L. Greene, S. 

G. Karabashev, D. A. Shulyatev, A. A. Arsenov and Y. 

Mukovskii: Phys. Rev. B 56 (1997) 13705. 
[13] L. Vasiliu-Doloc, J. W. Lynn, Y. M. Mukovskii, A. A. Arsenov 

and D. A. Shulyatev: J. Appl. Phys. 83 (1998) 7342. 
[14] C. A. Ramos, M. T. Causa, M. Tovar, X. Obradors and S. 

Pinol: J. Mag. Mag. Mater. 177-181 (1998) 867. 
[15] S. E. Lofiand, V. Ray P. H. Kim, S. M. Bhagat, M. A. Man- 

heimer and S. D. Tyagi: Phys. Rev. B 55 (1997) 2749. 
[16] C. V. Mohan, M. Seeger, H. Kronm[u]ller, P. Murugaraj and 

J. Maier: J. Mag. Mag. Mater. 183 (1998) 348. 
[17] A. Schwartz, M. ScheflJer and S. M. Anlage: Phys. Rev. B 61 

(2000) 870. Note that this experiment is for x = 0.2. 
[18] K. Ghosh, C. J. Lobb, R. L. Greene, S. G. Karabashev, D. A. 

Shulyatev, A. A. Arsenov and Y. Mukovskii: Phys. Rev. Lett 

81 (1998) 4740. 



